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Abstract 

The operator content of the Baxter- Wu model with general toroidal 
boundary conditions is calculated analytically and numerically. These 
calculations were done by relating the partition function of the model 
with the generating function of a site-colouring problem in a hexagonal 
lattice. Extending the original Bethe-ansatz solution of the related 
colouring problem we are able to calculate the eigenspectra of both 
models by solving the associated Bethe-ansatz equations. We have 
also calculated, by exploring the conformal invariance at the critical 
point, the mass ratios of the underlying massive theory governing the 
Baxter- Wu model in the vicinity of its critical point. 

1 Introduction 

The Baxter- Wu model is the simplest non-trivial spin model with three- 
interactions. Its Hamiltonian is defined on a triangular lattice by 

H = — J OiO,jO k , 

<ijk> 
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where the sum extends over the elementary triangles, J is the coupling con- 
stant and Oi = ±1 are Ising variables located at the sites. Historically this 
model was introduced by Wood and Griffiths in 1972 [|TJ, as an example of a 
model exhibiting an order-disorder phase transition and not having a global 
up-down spin reversal symmetry. This model is self-dual |l], with the same 
critical temperature as that of the Ising model on a square lattice. In 1973 
Baxter and Wu || related the partition function of this model, in the ther- 
modynamic limit, with the generating function of a site colouring problem on 
an hexagonal lattice. Solving this colouring problem through a generalized 
Bethe-ansatz they calculated the leading exponents H] a = 2/3, fi = 2/3 
and 7] = 1/4 of the Baxter- Wu model. The equality of these exponents with 
those of the 4-state Potts model [4-6] added to the fact that both models 
have the same four-fold degeneracy of the ground state, induce the conjec- 
ture that they share the same universality class of critical behaviour. This 
conjecture is interesting since contrary to the Baxter- Wu model the 4-state 
Potts model is exactly integrable only at its critical point. However from nu- 
merical studies of these models on a finite lattice it is well known that both 
models show different corrections to finite-size scaling (at the critical point). 
Whereas in the 4-state Potts models [7-10] these corrections are governed by 
a marginal operator, producing logarithmic corrections with the system size, 
which brings enormous difficulty to extract reliable results from finite-size 
studies of the model, this is not the case in the Baxter- Wu model [11-13]. 

Nowadays with the developments of conformal invariance applied to criti- 
cal phenomena [13] the classification of different universality classes of critical 
behaviour becomes clear. Two models belong to the same universality class 
only if they share the same operator content, not only the leading critical 
exponents. The operator content of the 4-state Potts model was already 
conjectured from finite-size studies in its Hamiltonian formulation [|nj] and 
can be obtained by a Z(2) orbifold of the Gaussian model in a special com- 
pactification radius (see []n| for a review). In this paper we numerically and 
analytically calculate the operator content of the Baxter- Wu model with sev- 
eral boundary conditions. Part of our numerical calculation was announced 
fL2| . In order to do this calculation we generalize the original Bethe- 
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ansatz solution of the related site-colouring problem with periodic boundary 
condition. This is necessary because this solution only gives part of the 
eigenspectrum of the associated transfer matrix. We also extend our Bethe- 
ansatz solution in other cases where the site-colouring problem is not on a 
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periodic lattice. This extension enables us to obtain the operator content 
of this model and the Baxter- Wu model for more general toroidal bound- 
ary conditions. Our numerical study was done by numerically solving the 
Bethe-ansatz equations and the analytical work was done by studying these 



equations using standard techniques based on the Wiener-Hopf method [16 



We believe that beyond the Ising model the Baxter- Wu model is the 
simplest spin model that can be solved exactly for arbitrary temperatures. 
We explore this solution in order to obtain the mass spectrum of the massive 
theory describing fluctuations near the critical point. 

The layout of this paper is as follows. In section 2 we introduce a site- 
colouring problem on an hexagonal lattice whose generating function is ex- 
actly related, in the bulk limit, with the Baxter- Wu model. Our construc- 
tion is valid for some toroidal boundary conditions, generalizing the original 
construction [tl for periodic lattices. In section 3 we present the transfer 
matrix associated to the colouring problem and calculate its eigenvalues by 
the Bethe-ansatz. In section 4 the operator content of the Baxter- Wu model 
and the related site-colouring problem is obtained. The mass spectrum of 
the massive field theory governing the thermal and magnetic perturbations 
is calculated in section 5. Finally our conclusions are present in section 6 
and the analytical calculation of some of the conformal dimensions of the 
site-colouring problem is presented in an appendix. 



2 The Baxter- Wu model and the related site- 
colouring problem 

In this section we relate the Baxter- Wu model with general toroidal boundary 
conditions with a site- colouring problem (SCP) on the hexagonal lattice. The 
construction presented here generalizes those presented by Baxter and Wu 
[||] for the periodic case. 

Let us consider a triangular lattice with L (M) rows (columns) along the 
horizontal (vertical) direction, respectively. For convenience we take L as 
a multiple of 3, and decompose the lattice in three triangular sublattices 
formed by the points denoted by O, □ and A in Fig. 1. We attach at each 
site of a sublattice Ising variables {cr°}, {a n } and {<r A }, and the Baxter-Wu 
model, with the Hamiltonian (|1|) is given in terms of the simplest three-body 
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interactions with one spin in each sublattice. 

The model has a non-local Z(2) x Z(2) symmetry corresponding to the 
global change of variables in two of the sublattices, i.e., 

(T A -> aa A a° -> 6a° a D -> afea n (a = 6 = ±1). 

This fact implies that at sufficiently low temperatures the model have a four- 
fold degenerate ground state. Creating, at low temperatures, domain walls 
connecting those different ground states we see that in fact the symmetry of 
those long-range excitations is -D(4), like in the 4-state Potts model. These 
reasonings suggest that critical fluctuations of both models are described in 
terms of the same quantum field theory. In order to do a detailed study 
of several sectors of this underlying field theory we consider the Baxter- 
Wu model with general toroidal boundary conditions compatible with its 
Z{2) x Z{2) symmetry: 

= a °trL,j = bcr° +Ltj a°j = abaf +LJ , (2) 

where a,b — ±1. In Fig.l the simbols O, □ and A denote the bordering sites 
which are related by (2) to the bulk ones (filled simbols). Imposing periodic 
boundary condition along the vertical direction the partition function of the 
model Zf^ M can be written as 

Z^ M = Tr(f™) M , (3) 

where is the associated row-to-row transfer matrix. Its elements 

(oi, <Tjv \TP ^ \a[, a' N ) are given by the Boltzmann weights generated by 
the spin configurations {ax, <tat} and {a[, a' N } of adjacent rows and are 
given by 

(<7i, (TArlfgjK, a' N ) = exp U{J2 a j a 3+^( a j+^ + a 'j+l) | > ( 4 ) 

with K = -j^j, and on the right side the appropriate boundary condition 
(a, b) is taken into account. 

Following Baxter and Wu we can relate Zf^ M to the partition function 
or generating function Zff^ p M of a SCP on an honeycomb lattice with N = ^ L 



4 



rows and M columns. In order to do this, it is convenient to define link 
variables {A} at the links of the hexagonal lattice formed by the sublattices 
O and □ (heavy lines in Fig.l). These variables are given by the product of 
the site variables at the ends of the link (A = a D a°). In terms of these new 
variables the Hamiltonian is given by 

H=-J y £of(\ 1 +\ 2 ..- + \ 6 ), (5) 

ieA 

where the summation is only over the sublattice A and Ai, A2, \ are the 
link variables surrounding a given site variable cr A . Taking into account 
the boundary conditions (|2|) the partition function, in terms of these new 
variables, can be written as 

Zlxm = £ n{ ex p{^ A (/i A i + /2A2 + --- + /6A 6 )}^(l + A 1 A 2 ...A 6 )j, 

(6) 

where the product extends over all the elementary hexagons formed by the 
sites in sublattices □ and O, and surrounding a site variable <r A . The factor 
|(l+AiA2-..A6) in (^j) is necessary since the variables {A} are not independent. 
The factors fi (i = 1, 2, 6) in @ are constants defined on the links of the 
hexagonal lattice (Q-O) and depend on their relative location (see Fig. 1). 
If both ends of a link i (with corresponding variable Aj) belongs to the bulk 
of the lattice (•-■) fi = 1, if one of its ends is on the border (O-m or • — □) 
fi = ab, and if both ends are on the border (O-D) fi = a. To proceed it is 
convenient to rewrite the last product in (|), surrounding each lattice site 
A as 

i(l + A 1 A 2 ...A 6 )= £ g(Xi, fi A )g(\ 2 , fi A ) ■ ■ ■ g(\ 6 , fi A ), (7) 

where 

^(A,^) = 2- 7 / 6 (A + /i+|A- yU |). (8) 

Substituting the expression (^) into equation @ the partition function will 
be expressed in terms of a single sublattice A with site variables {<J A } and 
{/i A }, and link variables {A}. Taking into account the boundary factors it is 
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straightforward to write 



{a A ,fi A } <i,j> {\i,j=±l 

(9) 

where < i,j > are links on the hexagonal lattice formed by the sublattices □ 
and O (see Fig. 1). In order to derive @ we are forced to restrict ourselves 
only to boundary conditions (|2|) where a = b = ±1. Fortunately this does 
not restrict our analysis since due to the D(4) symmetry of the model, the 
eigenspectra of the transfer matrices with boundary condition a / i and 
a = b = — 1 are degenerate. Summing over {A} in (§) we obtain 

z l™m= Yl II w ( ff f.ft; ff f^f). ( 10 ) 

{a A ^ A } <i,j> 

where 

w(cr h Hi\ VjiHj) = 2 ~ 1/3 [exp (Kfa + <7j)) + fiiHj ex P (r K (?i + a i))\ ■ ( n ) 

Finally, following Baxter and Wu || we now associate the above partition 
function to the generating function of a SCP with 8 colours. We associate 
odd colours 1, 3, 5 and 7 at a given site A according to the values of the 
variables (a A ,/i A ) attached at the site: 

(++)-! ("+)-3 (— )-5 (++)-7. (12) 

The relation (O) tell us that links connecting colours 3,7 and 1,5 should 
be forbidden in the SCP. This constraint can be easily implemented J| 
by introducing the even colours 2, 4, 6 and 8 on the sublattice O, forming 
with the sublattice A an hexagonal lattice with N = |L (M) sites in the 
horizontal (vertical) direction. The constraint is that all nearest-neighbour 
colours on this hexagonal lattice must differ by ±1 (modulo 8). For M — > oo, 
the partition function Zf^ is then related to the generating function Z^m, 
given by 

y SCP _ m n 2 ~"8 _ yBW 

^N,M — Z l Z 2 ' ' ' z 8 — ^ L,Mi \ LO ) 
{C} 

where n\,U2, ■■■,ns is the number of sites coloured with colour 1,2,.. .,8 in a 
given configuration C. The fugacities Zi (i = 1,2,. ..,8) are obtained from 
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(ED 

zi = z 3 = z 5 = z 7 = 2sinh(4K), , . 

z^ 1 = z 4 = Zq 1 = z 8 = sinh(2i^) = t. 

The critical point of the Baxter- Wu model and of the SCP is given by the self- 
dual point t — t c — 1 [Q. Concerning the boundary conditions the relations 
are as follows. The with periodic boundary condition (a = b = 1 in 

(0)) is related with Zff^, also with periodic lattice. The boundary condition 
a = b = — 1 in (2) is related to a boundary condition in the SCP such that 
if at site we have a colour variable cij = 1, 3, 5 or 7 at site (iV + l,j) 

the colour variable should be c/v+ij = 3, 1, 7 and 5 respectively. We do not 
need to mention even colours since N is even (see Fig. 1). 

3 The transfer matrix and the Bethe-ansatz 
of the SCP 

In this section we derive the row-to-row transfer matrix of the SCP and 
generalize its Bethe-ansatz solution presented by Baxter and Wu 0. The 
solution presented in ||, as we shall see, only gives part of the eigenspectra 
of the transfer matrix of the SCP, with periodic boundary conditions. Here 
we extend the solutions for the periodic case and also derive the Bethe-ansatz 
equations for the more general boundary condition 

= Ci + N,j + 2k (mod 8); k = 0, 1, 2 or 3, (15) 

along the horizontal direction. The periodic case corresponds to k = 0. 
Following |J, for a given configuration {cjj} of colours on the hexagonal 
lattice, we say we have a dislocation on a given link of the lattice wherever 
the colour on the right end of the link is smaller (modulo 8) than that on 
the left end. In Fig. 2 we show two configurations with the corresponding 
dislocations (dotted lines) for a lattice with width iV = 4. In Fig. 2a the 
lattice is periodic (k = 0) and in Fig. 2b the boundary condition is given by 
(13) with k — 1. 

The generating function of the SCP can be written as 

Z S ^ = Tr(f^) M , (16) 
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where T^ p is the associated row-to-row transfer matrix. This transfer 

matrix has elements ({C}\T^ P \{C}) given by the product of the Boltz- 
mann weights due to the colour configurations {C} = {ci, C2, .., cn} and 
{C} = {c[, c' 2 , d N } of two adjacent rows. If the configuration produced 
by {C} and {C} contains only colours differing by ±1 (modulo 8), we have 

N 

({c}\fffi p \{C}) = n(*(*)*W)) 1/a , (17) 

i=l 

with fugacities Z{ defined in flT4|). On the other hand if the configuration does 
not satisfy this constraint, ({C}\T^ P \{C}) = 0. It is simple to see that 
this requirement implies, for arbitrary values of k in (|i"5|) , the conservation of 
the number n of dislocations along the vertical direction. Consequently the 
Hilbert space associated to T^ p can be separated into block- disjoint sectors 
labelled by the values of n. The possible values of n depend on the boundary 
condition fll5|) . For k = 0, 2 (1, 3) they are even (odd), and are given by rij, 
where 

< U j = N- K-ij < 2N j = 0, ±1 ± 2, ... (18) 

The colour configurations {C} and {C}, in a sector with n disloca- 
tions, can be conveniently expressed by the sets (m;X) = (m; Xi,X2, x n ) 
and (m'; X') = ( ), respectively. The odd numbers m and 

m! give the colour at the first site and the sets X = (x\, X2, x n ) and 
X' = ( x'i, x' 2 i x' n ) give the position of the dislocations on the row. The sets 
X and X' should satisfy 

1 < x x < x 2 < ■ ■ ■ < x n < N , 1 < x[ < x' 2 < ■ ■ ■ < x' n < N, (19) 

and should have no more than one repeated value of x or x' if they are odd, 
or more than three repeated values if they are even. In Fig. 2a we show 
the configurations (m;X) = (1;2,3,4,4) and (m';X') = (3; 2, 2, 3, 4) which 
belong to the sector with n = 4 and periodic boundary condition. In Fig 
2b we show the configurations (m; X) = (1; 2, 3, 4) and (m'; X') = (3; 2, 2, 2) 
belonging to the sector with n = 3 in the lattice with boundary condition 



K — 1 in (p~5[ ) and width N = 4. Two configurations X = (xi,X2, ■■■ } x n ) and 
X' = ( Xi, x' 2 i -i x' n ) are connected through the operator T^ p if beyond (|19]) 

x'j = xj - 1, if Xj odd ^ 2Q ^ 

00 OC j ^ tJC rj 1 > 00 if 00 j G'VGH • 
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In the case Xi = Xi+i we have an additional constraint x\ 7^ x' i+l . We should 
notice, for arbitrary boundary condition, the identification 

X = (0, x 2 , x 3 , ...,x n ) = (x 2 , x 3 , x n , N). (21) 

The transfer matrix, in a given sector with n dislocations is now given by 

i-(-i) K n 

(m,X\f^ P \m',X'} = z m+ { Yi w ( m + X J + x j ~ 2 j), w ( m ) = ( z mZ m+ i) 1/2 , 

3=1 

(22) 

if X and X' satisfy (|1^) and (pO]), and is zero otherwise. In the sector with n 



dislocations the eigenvectors ip^ of T^ p , with eigenvalue A can be written 

as 

^ {n) = V1 8 fm(x u x 2r ■ ;x n )\m;X), (23) 

{m,X} 

where the summation is restricted to the configurations with n dislocations, 
and f m (xi, x 2 , ■ ■ •, x n ) are unknown amplitudes. The eigenvalue equation for 
Tfik F is given by 

* / n \ 

£ I] w ( m + X J + 4 - 2i) / m+2 (X') = A/ m (X), (24) 
x' \j=i J 

where the asterisk indicates that X and X' satisfy the conditions fli~9|) and 
fl2"0p. The relation ( ]2"T| ) implies that in ( ^ ) the amplitudes having x[ = 
should be replaced by the boundary condition 

f f f i-(-i) K f / / 

f m+2(0, x 2 , £3, • • •, x n ) = z m+ i f m (x 2 , x 3 , ■ ■ •, x n , N). (25) 



Due to the values of the fugacities (JT^) it is simple to see that T^k , 



besides conserving the number of dislocations, also has an additional Z{2) 
symmetry (eigenvalues e = ±1), since adding 4 (modulo 8) to all colours in 
a given configuration does not change its weight in the generating function, 
that is 

f m+ 4(X) = ef m (X). (26) 

Following Baxter and Wu U we assume the following Bethe-ansatz for 
the amplitudes 

f m (X) = J2a(P)(l>Pi(rn-2,x 1 )---(f>p n (m-2n,x n ), (27) 
p 
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where the summation is over all the n! permutations P = {Pi, P2, P n } of 
integers {1,2, ...,n}. We require the existence of n wave numbers kj (j = 
1, 2, n) and signs €j (j = 1, 2, n) such that 

^(m,x) = e^(m + 4, re) = { a ^^kx) x odd 
J J JV I Oj^exp^i/CjX), x even. 

Observe that the Z(2)-parity eigenvalue of the wave function is given by 
e = ]lj=i e j) an d it is even or odd depending on the numbers of negative 
values of €j. The Bethe-ansatz solution presented by Baxter and Wu || only 
gives the symmetric eigenvalues (e = 1), for periodic boundary conditions 
(k = 0). We can follow the same procedure as in || in order to derive the 
Bethe-ansatz equations. We have to consider various possible choices of X 
to determine the eigenvalue A. 

Firstly let us consider the case where all dislocations are located at dis- 
tinct positions, i.e., x± 7^ x 2 7^ • • • 7^ x n . Eq. ([2~4|) is then replaced by 

'^2w(m + x + x')(f>j(m + 2,x') = Xj(j)j(m,x), j — l,...,n, (29) 

x> 

where we have denoted 

A = Ai • • • A n . (30) 

Actually Eq. (|29|) represents two equations corresponding to x odd or even 
and can be written as 

^j,m+2^ / j,m+2 = ^jVj,mi (31) 

where 



\ '~ x j,m— 1 ■ rl -j,m / \ "jjm / 



with 

ctj,m = w(to) exp(— ifcj); Aj >m = w;(m) + w{m — 2) exp(— 2ikj). (33) 
Using the fact that 

( a j,m-2 J _ £ . f a j,m+2 \ 
\ i>j,m-2 ) 3 \ bj,m+2 J ' 

it is simple to see that Xj can be obtained from the eigenvalue equation 

(Tj tm Tj^ m+ 2)Vj,m+2 = X^ejVj jTn+ 2- 
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Solving for the eigenvalues of Tj jm Tj jm+ 2 we see that Xj does not depend on 
the value of m (as it is expected) and satisfies 

A* - e j A J 2 [exp(-z4fc i ) + A exp(-2zA; i ) + 1] + exp(-4zA; i ) = 0, (34) 

where A = t+l/t. The solution of fl34|) is given by Xj = ^/e] exp(ej — ikj ), 
with 

ef ] = 1/2 In (xj + sjyjx? - l) , Xj = cos(2A;f j) ) + A, (35) 
and Sj = ±1. Equation (|31~1) gives the relations 

where 

%,m = ejyftjAj^Xj - eja) m _ x l Xj 1 ). (37) 

It is important to verify that f2j im+ 2fij, m = 1, so that at t — t c — 1, = 1. 

Secondly let us consider the case where two even dislocation positions x 
coincide. For convenience suppose X\ = x 2 = x (even). If we require that the 
ansatz (p7"D satisfies (|4]) with eigenvalue given by (|30D and ( p5|) the equation 

^ iu(m + x + x' x - 2)w(m + x + x' 2 - 4)/ m+2 (x' 1 , x 2 ) = AiA 2 / m (x, x), (38) 



x l > a '2 



must be fulfilled, where 

f m +2(x 1 ,x 2 ) = a(l,2)0i(m,xi)0 2 (^-2,x 2 ) + a(2, l)0 2 (m, Xi)0i(m - 2, x 2 ), 

(39) 

and the asterisk in Eq. (|31f ) indicates summation over the possible configu- 
rations (x[, x' 2 ) = (x — 2, x — 1), (x — 2, x), (x — 1, x). Using (|28|) , (|32|) , (0) 
and (|36| ) in ( ffSj ) we obtain after some algebra (the same algebra as in Q), 
that the ratio 

S ^2) _y^cosh(e^+»fc 2 ) = VH 

a(2,l) v^cosh^+iJfeO v 7 ^ 1 ^ 

is also independent of m 
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More generally in order for the Bethe-ansatz fl2"T|) to work we should have 

~ = a{...,j,l,...) = y/ej cosh(e^ ) + ikj) = 
jl ~ a(...,l,j,...) v ^cosh( e ; si) +2%) v^7 jh 

for all permutations of adjacent elements in a(P). Following the same steps 
as in || it can be proved that equations ([IT]) are enough to ensure the effec- 
tiveness of the ansatz ( p7|) in the case of triple coincidence of dislocations. 
In order to complete the solution we still need to fix the wave numbers 
{ki, k 2 , k n }. As usual this is done by inserting (|27|) in the boundary con- 



dition (gSJ), 

i-(-i) B 



£a(Pi, ■■■,P n )<t>Pi(m,0) [<j)p 2 {m - 2, x 2 ) ■ ■ ■ </>p n (m - 2n + 2,x n )] = z m+ \ 
J2a(P 2 ,...,P n ,Pi) [4>ih{rn - 2, x 2 ) ■ ■ ■ (j> Pn (m - 2n + 2, x n )\ Pl (m - 2n, N), 

(42) 

where in the right-hand side a circular shift on P was done. This equation is 
fulfilled if we require that 

i-(-i)" 



a(P l ,P 2 ,...,P n )4> Pl (m,0) = z m+ f a(P 2 ,P 3 ,...,P n ,P 1 )0 Pl (m-2n,N). (43) 
Using d|), d|) and © we then obtain 

n 

exp (iNkj) = -Q Bjl > j = 1, 2, n (44) 



where 



tat^XX-i] and (%=(P*±) 4 ■ (45) 



Actually we have in (|44j) two distinct sets of equations, a first one for m = 1 
and m = 5 and a second one for m = 3 and m = 7. These sets must be solved 



simultaneously for the validity of the Bethe-ansatz (P7|). From Eqs. (14) and 
(B3)-(37) we verify that those equations degenerate (p^\ = 1, m = 1,3, 5,7; j = 



1, 2, n) in the following cases: (a) k = or 2 for arbitrary values of temperatures 
i, (b) k = 1 or 3 only at the critical temperature t = t c = 1. In both cases the 
Bethe-ansatz equations are given by 

n 

eM-iNkj) = -(-Ff-Vin B H , e = ±1 (46) 

l=i 
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with Bj t i and n_ given by (41) and respectively. For a given value of e = ±1 
the prefactor 5- = (— l) n ~, in Eq. (f4q), may be positive or negative depending 
on the particular choice of the set {ei,...,e n }. If n is even these two choices 
are equivalent. The solution of both equations are the same except that one of 
the quasi- moments differs by the value ir (mod 2ir). However if n is odd the 
situation is different, and we should consider both equations. Their solution gives 
us independent wave-functions. 



4 The operator content of the Baxter- Wu 
model and the SCP 

In this section, by exploiting the conformal invariance at the critical point, we 
are going to derive the operator content of the Baxter- Wu model and the SCP. 
The conformal anomaly c and anomalous dimensions (x%,X2, • ••) are obtained in a 
standard way from the finite-size behaviour of the eigenspectra of the associated 
transfer matrix, at the critical temperature. If we write T = exp(—H), then in 
a strip of width L with periodic boundary conditions the ground-state energy, 



Eq(L), of H behaves for large L as fl7[| 

E (L) 



+ o(L- 2 ), (47) 



L °° 6L 2 

where is the ground-state energy, per site, in the bulk limit. Moreover, for each 
operator O a with dimension x a there exists a tower of states in the spectrum of 
H with eigenenergies given by [13, ^] 



E^ m ,(L) =E + ^(x a + m + m') + o^" 1 ), (48) 

where m,m' = 0,1,2,... . The factor v s appearing in the last two equations 
is the sound velocity and has unit value for isotropic square lattices. We can 
calculate directly the higher eigenvalues of and T?J? P of our model by a 

numerical diagonalization, using for example the power method. However since 
these matrices are not sparse and have dimension 2 L , for a horizontal width L, 
we cannot compute the eigenspectra by numerical diagonalization methods for 
L > Lq ~ 18 in the case of the Baxter- Wu model. 

Instead of a direct calculation we can explore the relation (|l^) among the 
Baxter- Wu model and the SCP, and solve numerically the associated Bethe ansatz 
equations derived in section 3. If we write T = exp(-H) for both models, the 
relation (13) implies 

Tr{e- MH L W ) = Tr(e- MH N CP ), (49) 
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where N = Its is important to observe that although H PW and Hfp p have 
the same dimension 2 L they may have different eigenvalues. Indeed that is the 
case, specially if t 7^ t c = 1 as we can verify by a brute force diagonalization of 
these transfer matrices on small lattices. In order to present our results let us 
rewrite the Bethe-ansatz equations (^) in a convenient form. We are going to 
choose the prefactor 5- = (— l) n ~ = 1, since as we discussed in the last section, if 
n is even all the energies can be obtained from a given choice of <5_ , and for n odd 
the energies obtained by different choices of 5- are complex-conjugate pairs. The 
eigenvalues of Hff' p in the sector with n dislocations are given by 

4«> = -^ln(16t 2 (l +i 2 )) - f>5 Sj) " ikf j) ), (50) 

i=l 

where 



e ( f j) = 1/2 In (xj + s^x) - l) , Xj = cos(2k ( f j) ) + t + 1/t, (51) 
with 1 = si = S2 = ■ ■ ■ = s n ^i = — s n _i + i = . . . = —s n , and I = 0,1, ...,n fixed. 

(s) 

The quasi-momenta {k- 3 } are obtained by solving the equations 



exp 



p =i ycosh(ep ' +ikj 1 1 



where e = 1 (e = — 1) gives the even (odd) part of the eigenspectrum, with respect 
to the Z(2) symmetry of the SCP discussed in the last section. Strictly speaking 
these energies may only give part of the eigenspectra, since the completeness of the 
Bethe-ansatz solution presented in the last section is an open question. Numeri- 
cally we have studied these equations extensively for lattice sizes up to N ~ 200 and 
part of our results at t = t c = 1 were presented in [12]] . For example, the ground- 
state energy for Hf j CP corresponding to the boundary conditions k = 0, 1, 2 and 3 



given in (15) belongs to the sector where n = N — k, s\ = S2 = ■ ■ ■ = s n = 1 and 
e = 1. In table 1 we present these ground-state energies, per site, for k = 0, 1 and 
2 (k = 3 is degenerate with k = 1). The conformal anomaly is obtained by using 
([47|). The bulk energy e^ p = — |ln6 can be obtained from the exact solution 
in the bulk limit j|] and the sound velocity Vg = \/3/3, can be inferred from 
(flg|) and an overall analysis of the dimensions appearing in the model. With these 
values the ground-state energy (first column in table 1) gives us the estimators 
c(N) presented in table 2. As expected the conformal anomaly is c = 1, as for the 
4-state Potts model. The direct calculation of the eigenspectra of H SCP and H BW 
for small chains shows us that although the eigenspectra of both models is not the 
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same, several eigenvalues coincide. This is the case for the ground-state. Conse- 
quently by using the bulk limit value || = — ^ln6, and the sound velocity 
vf w = y/S/2, we obtain the expected value c = 1 for the Baxter- Wu model. 

The dimensions defining the operator content of the model are obtained from 
the large-N behaviour of the energies of excited states. Let us concentrate on the 
SCP. The Bethe-ansatz equations fl52] ) are the same for all boundary conditions, 
specified by k (0,1,2 and 3) in (|l~5|), only the allowed values of n, given by (|ilj|) 
depend on the particular value of k. Using (fl8|), the finite-size sequences for some 
dimensions are shown in table 3. In this table x e j(N — n,l) are the estimators of 
the dimensions associated to the jth lowest eigenenergy in the eigensector labelled 
by n and e. The values of I used in (|50|)-(|52|) to obtain the corresponding energies 
are also shown. The numerical solution of (|50|)-(|52|) was done by the Newton type 
method. The roots {kj} for all the solutions that we obtained are real. Although 
we cannot discard this possibility for higher states, we do not find any string- 
type solution. The numerical solution of (|50|)-(|5"2"|) is not easy for N ~ 200 due 
to numerical instabilities, and some tricks are necessary. In most cases we solve 
initially these equations for small values of t, where a good guess can be given, and 
use the solution obtained as the initial guess for a larger value of t. We repeat this 

process up to t = t c = 1. As a result of our extensive calculation of the eigenspectra 
oiH SCP 

we arrive at the following conjecture. Namely, the dimensions of primary 
operators in a given sector labelled by n = N + k + 4 j (j = 0, ±1, ±2, ...) of the 
H SCP with boundary condition k (0,1,2 or 3) are given by 

^, 9 = ^(V + j), P = 3-j, 9 = 0,±1,±2,... . (53) 



In obtaining ( |53| ) , the value of Eq in fl48|) is the value of the ground-state energy for 
the periodic (k = 0) SCP. Moreover the number of descendants with dimensions 
Xp^q +m + m' (m, m G Z) is given by the product of two independent Kac- Moody 
characters. This allows, by using (^7]) and (|4"8|), to write the generating func- 
tion of the SCP with boundary condition k, up to order exp(— M/N) 
(M, N oo) as 

Zlf/ M ^) = eM-ei CP MN)z^Q\z) £ £ z^ + ^) (54 ) 

P £(Z-k/4) qeZ 



where 

z = exp(^f^) @(z) = (1 " a")" 1 - (55) 

n=0 

Since the Bethe-ansatz roots {kj} of ( |52|) are real numbers, we can apply analytical 
methods [21-24] based on the Wiener-Hopf method to obtain the finite-size cor- 
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rections of the eigenenergies. We calculate the finite-size corrections of the lowest 
eigenenergies in the sector with n dislocations and parity e. Since these calcula- 
tions are rather technical we present them in appendix A for the interested reader. 
These analytical results are in agreement with the conjecture (|53|), obtained from 
the numerical solutions of (|50|)-(|52|). 

Let us return to the Baxter- Wu model. Consider initially the model with pe- 
riodic boundary conditions (a = 6 = 1 in @). Comparing the eigenspectra of 
fflf and T£f? p , obtained by a direct diagonalization on small lattices, we verify 

that many of the dimensions Xp) q appearing in (|53| ) are absent. For example, the 
energies producing the estimator x^~(0, 1) in the 5th row of table 3 only appear 
in . Following for large lattices the energies which are exactly related in 

both models we verified that the lower dimensions in the Baxter- Wu model, with 
periodic boundary condition, are given by x = 0, |, |, |, . . ., and appear with de- 
generacy d x = 1, 3, 1, 9, respectively. Due to its D(A) symmetry the Baxter- Wu 
model has the same eigenspectra for the non-periodic boundary conditions given in 
(||), i.e., (a, b) / (+, +). We have shown, at the end of section 3, that the partition 
function in these cases is exactly related with a SCP with boundary conditions 
not included in fll5|). Actually the application of the Bethe-ansatz in this case, 
if possible, is more difficult since the number of colours in a row is not a good 
quantum number anymore. However at t = t c = 1 our direct calculations of the 
eigenspectra on small chains shows that there exist exact coincidences between the 
eigenvalues of TP \ and those of T?9 P , which are given by the Bethe-ansatz equa- 
tions (|5C|)-(|5"2j). These coincidences enable us to verify that the lower dimensions 
of the Baxter- Wu model with boundary condition (a, b) ^ (+,+) in @ are given 
by x = 1/8,1/2,10/16,... and appear with degeneracy d x = 1,1,4,... These are 
the same dimensions reported in [14] for the 4-state Potts model with antiperiodic 
boundary condition. These results supplemented with the global eigenspectrum 
calculated for small systems, indicate that the operator content of the Baxter- Wu 
model is the same as that of the 4-state Potts model [|f|] and is given in terms of 
a Z(2) orbifold [|H|] of the Gaussian model. 

Before closing this section, since we have calculated the eigenspectra of H BW 
and H SCP for large lattices we can also calculate the dimensions of the operators 
responsible for the corrections to finite-size scaling in both models. Since these 
calculations were already presented earlier (see Eqs.(ll)-(12) and Table 3 in Jl^]), 
we only mention that x 7 = 4 is the lowest dimension of the operator responsible 
for the finite-size deviations of the critical behaviour. This means that relations 
(p?^) and ( p8| ) have corrections which are power-like with the system size L. These 
corrections are like those of the Ising model and different from those of the 4- 
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state Potts model. This explains why the finite-size studies of the Baxter- Wu 
model have good convergence, in contrast with the 4-state Potts model, where 
the operator responsible for these corrections is marginal (x 7 = 2) producing 
logarithmic behaviour with the system size. 

5 The off-critical properties of the Baxter- 
Wu model 

The Baxter- Wu model and the SCP have a massive spectra at t ^ t c = 1. A 
continuum field theory describing the long-distance physics in this phase can be 

obtained in the neighborhood of the perturbing thermal parameter 5 = t — t c ~ 0. 

Such a field theory will be massive and the masses can be estimated from the 
finite-size behaviour of the eigenspectra of H = — In T. We can calculate the mass 



spectrum by applying the scheme proposed by Sagdeev and Zamolodchikov [23] 
in the study of the Ising model under the influence of magnetic perturbations. 
According to this scheme we should initially calculate the finite-size corrections 
of the zero-momenta eigenenergies E k (5,L), k = 0,1,2,..., at the conformal in- 
variant critical point 5 = 0. From our analysis presented in the last section these 
corrections are governed mainly by an irrelevant operator with dimension x 7 = 4 
and have integer power-law behaviour with the system size L. According to con- 
formal invariance |l8|] E k (5,L) should behave as 



E k (L) = eoo L + [x k - — ]+ a k ^L~ 3 + a k)2 L~ 5 + ■■■, (56) 



L V 12 



where x k is the conformal dimension associated to E k and a k> % (i = 1,2,...) are 
L-independent factors. According to the scheme of [23|, if the perturbed operator 



which produces the massive behaviour has dimension y, we should calculate the 
eigenspectra in the asymptotic regime 5 — > 0, L — > oo, with 



i 



X = 5~L (57) 

kept fixed. In this regime (|5^) is replaced by 

E k (6,L) = e 00 L + 6^F k (X) + a k>1 5^G k {X) + a k>2 6^H k (X) + --- . (58) 

The masses of the continuum field theory are obtained from the large-X behaviour 
of the functions [^] F k (X), and are given by 

m k ~ F k (X)-F (X), (59) 
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where Fq(X) is associated in ( pq ) with the ground-state energy. 

In the present application, the thermal fluctuations are produced by the energy 
operator, which has dimension y = x e = 1/2. Since we are going to calculate the 
eigenenergies of the Baxter- Wu model by exploiting its connection with those of 
the SCP, it is important to compare their eigenspectra for small lattice sizes. In 
table 4 we represent for t < t c (T > T c ) the relative location in the eigenspectra 
of the lower zero-momentum energies of both models. The eigenenergies in the 
same line are exactly degenerate on the finite lattice. In this table we also show 
the parity-quantum number e ± 1, of the eigenenergies of H SCP . 

It is important to mention that although the Baxter- Wu model and the SCP 
are exactly related, the parameter t has quite a different effect in both models. 
In the case of the Baxter- Wu model it drives the system from an ordered phase 
(T < T c , t > t c = 1) to a disordered phase (T > T c , t < t c = 1). On the other 
hand, as we can see form (14), in the SCP it drives the model from an ordered 
phase rich in colours 4 and 8 (t > t c = 1) to another ordered phase rich in colours 
2 and 6 (t < t c = 1). This fact implies that even for T > T c we should have 
in the SCP an infinite set of states, including the ground state, that degenerate 
exponentially with the system size (E n — E' n ~ exp(— aL)), which certainly is not 
the case for the Baxter- Wu model in its disordered phase. In table 4 the pair of 
levels Ef and E-^[ p (i = 0, 2 and 4) degenerate exponentially. Baxter and Wu 
Q in their original calculation of the exponent a of the Baxter- Wu model used 
the excited energy E2 CP (e = 1), instead of E^ CP (e = —1) = E PW . However, as 
we mentioned, these energies become exponentially degenerate with system size, 
not changing their exact result a = 2/3. 

Exploring the correspondences presented in table 4 and using (|57|)-(|59|) we 
can calculate the mass ratios of the underlying massive field theory governing the 
Baxter- Wu model for T ^ T c . They are calculated from the asymptotic regime 
X — > oo of the finite-size sequences 

F k (X,L)-F (X,L) m k 
R ^ L ^~ Fl{X ,L)-F (X,L)^W (60) 

The functions F^X, L) are obtained by using in ( ^8| ) the finite-size sequences of 
the zero-momentum states (fc=0,l,2...). The exact degeneracy of E PW (see table 
4) implies the first mass mi is triple generated. From the equality E PW = Ef CP 
we can calculate the second mass m<i by solving the Bethe-ansatz equations derived 
in section 3 for the SCP . Unfortunately, although trying hard, we were not able 
to find the Bethe-ansatz roots that would correspond to this energy. However, 
applying the Lanczos method directly in H SCP we calculate this eigenenergy up 
to L = 21, in the Baxter- Wu model. In table 5 we show the estimators R2(X,L) 
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obtained by using in (Eq)-(E(j) L = 15, 18 and 21. These results are consistent with 
the conjecture m,2 = \/3mi- In the case of mass m 3 our results are more precise 
since we were able to calculate E§ w = E§ cp (e = 1) for lattice sizes up to L=150, 
by solving the Bethe-ansatz equations (|50|)-(^2|). In table 6 we show the estimator 
Rz(X,L) for some values of X, obtained by using in (|58|)-(|60|) L = 144, 147 and 
L=150. These results show clearly that 777,3 = 2m%. The numerical analysis of 
other higher energies in the spectrum shows that a continuum starts at 771,3. 

The mass ratios we obtained should be the same as those of the 4-state Potts 
model, since we expect both models share the same universality class of critical 
behaviour. In fact they coincide with the masses previously conjectured [^4|] for 
the 4-state Potts model and are also given by the masses of a sine-Gordon model 
1 25] at a special coupling, i. e., 



7T 

m i+ i = mism(-i), i = 1,2,3. (61) 

D 

To conclude this section we mention that we also studied the effect of mag- 
netic perturbations in the Baxter- Wu model. This was done by calculation of the 
eigenspectra of with the addition of an external magnetic field h. In this 

case 5 = h and y = 1/8 in (^), and the masses we obtained are consistent with 
those reported in [^] for the 4-state Potts model. However our results in this case, 
specially for larger masses, lack precision because we had to calculate directly the 
eigenspectrum of T^^i since the equivalence with the SCP presented in section 
2 is not valid anymore and unlike Tyv^M this matrix is not sparse. 



6 Conclusions and Comments 

The operator content of the Baxter- Wu model was calculated for several boundary 
conditions by exploiting the conformal invariance of the infinite system at the 
critical point. Our results are calculated analytically and numerically for very 
large lattice sizes. This was possible due to the relation between the Baxter- Wu 
model and the SCP. Actually we showed that the partition functions of both 
models are exactly related for several boundary conditions (see section 2) and 
we were able to extend the original Bethe-ansatz solution for most of these 
boundaries (see section 3) 

The operator content of the SCP with several toroidal boundary conditions 



(see (|53|)-(|55|)) is the same as those of a Gaussian model with dimensions [27] 



771^ 

x n ,m = gn 2 + — , (62) 

4y 
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where g = 1/8 is the compactification radius and n,m S Z are the vorticity 
and spin-wave number, respectively. However only part of the eigenspectra of 
both models coincide. Our analysis (section 4) shows that the dimensions of the 
Baxter- Wu model, for several boundary conditions are given by a Z(2) orbifold 
of the above Gaussian model. This operator content coincides with the 4-state 
Potts model, indicating that indeed both models share the same universality class 
of critical behaviour. It is interesting to remark that whereas for the SCP the 
operator content is given in terms of characters of the Kac-Moody algebra, in the 
Baxter- Wu model the characters are those of the Virasoro algebra. 

On the other hand a similar exact relation as that between the SCP and the 
Baxter- Wu model also exists between the 4-state Potts model and the 6-vertex 
model at its isotropic point (7 = 0), or equivalently the quantum XXX chain 
(anisotropy 7 = 0). The operator content of this model is given by a combination 
of the dimensions given in (|62| ) but with g = 1/2. In the 6-vertex model, or XXZ 
chain, the point where g = 1/8 corresponds to the anisotropy 7 = Sir/ 4, and is 
the so-called Kosterliz-Thouless point. This implies that exactly at the critical 
point (L — ► 00, T = T c ) the Baxter- Wu model and the 4-state Potts model are 
governed by the same conformal theory, but deviations from the critical point, like 
for example the finiteness of the lattice, will be governed by an effective Gaussian 
model with different compactification radius. It is known || that in the case of 
the 6-vertex model or XXZ chain the finite-size corrections are ruled mainly by the 
operator with dimension £0,2 in (p^ ) besides the descendant of identity operator 
with dimension 4. This implies the appearance of logarithmic corrections, with the 
system size, at g = 1/2 since xo,2 = 2 and the corresponding operator responsable 
for such corrections is marginal. On the other hand at g = 1/8 we only have integer 
power-law corrections with the system size, since in this case the operator with 
dimension 4 dominates the finite-size correction. This explains why although the 
Baxter- Wu model and the 4-state Potts model are given by the same Z(2) orbifold 
of a Gaussian theory, they show quite different behaviour at finite lattices. 

In section 5 we calculated the mass spectrum of the underlying field theory 
governing the Baxter- Wu model around its critical point. In the case of thermal 
perturbations we obtained the masses given in (|6l"D which are the same as those of 



the 4-state Potts model [24] and are also the masses of a special point of a massive 
sine-Gordon field theory p5| . Finally in the case of magnetic perturbation our 
numerical results, although poorer, are consistent with the same masses reported 



earlier for the 4-state Potts model |2£ 
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Appendix 



A Analytic calculation of the leading finite- 
size corrections 

In this appendix we calculate analytically the leading finite-size corrections for 
some of the eigenenergies of the SCP at the critical point t = t c = 1. We will 
calculate the finite-size corrections of the lowest energies E n {e) of the Hamiltonian, 
H SCP = —InTscp, in the sector with n dislocations and Z(2)-colour parity e 
(n = 0, 1, 2, ...; e = ±1). Since the associated roots of the Bethe-ansatz equations 
are real numbers our analytical calculations are based on the method pioneered by 
de Vega and Woynarovich [pi] and Hamer [20] and refined by Woynarovich and 



Eckle |21j (see also |22]). In order to apply this method it is convenient to change 
the variables {k\, k 2 , • k n } appearing in ( |4"6| ) into new variables {m, u 2 , u n } so 
that Bjj become a function of the difference Uj — u\. This was done by Baxter 
and Wu || for arbitrary temperatures and Bjj(iij — u{) are now given in terms 
of elliptic functions. At t = t c = 1 these elliptic functions become hyperbolic 
functions, with 

Bjj = — exp(— iO(uj — ui)) = — itanh(uz — Uj —m/4). (A.l) 



In terms of the variables u the quasimomenta k(u) and the factors e(u) in ( |5l|) are 
given by 

eW ^ 1/21n ( ^M±^| ). (A,) 
Vcosh(2nJ — v2/2 / 



The Bethe-ansatz equations (46) or (|>2|) can be written as 



I- 1 { 1 n 1 

_2 = _ ) k ( Uj) + _ ^ @ {u . _ Ul) K j = i, 2 , ...„, (A.4) 

where 21 j are integers or half-odd integers depending on the value of e. The values 
of Ij for the lowest eigenenergy in the sector with given values of n and e, which 
we are interested in, are 

t t t n—l — en — l — i n — l + e „ 1-e . . 

h,I 2 ,-,In = ^ , ^ + 1,-, ^ , e = ——. (A.5) 
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Following a standard procedure [JOI we define the density of roots 

*w = ^. < A - 6 > 

where 

^ («) = ^ {*(«) + ^£e(«-«j)}- (A.7) 

When JV — ► oo ( A,6| ) becomes an integral equation whose solution gives the bulk 
limit of the density of roots. In particular, in the sector n = N this density of 
roots is given by 

4 cosh(4-u) cosh(8n/3) 

o-oo («) = ,, a x . • (A.8) 

7r cosh(8uj + 1 

In this limit the energy per site is given by || 

e ^ = -jjln(6). (A.9) 

The difference between the energy per site and the density of roots and their 
bulk-limit values can be expressed by 

fpn poo 

-f ~ eif P = ~ J ^ f(v)S(v)dv, (A.10) 

and 

2vr 



<r n N (u)-a 00 {u) = -^ I p(u-v)S(v)dv, (A.ll) 



respectively, where 

1 n 

S(v) = -^S(v- Uj )-a^(v), (A.12) 
i=i 

3 smh(4u) 

„, \ , [°° sinh(7rx/8) ,. . , 

/(«) = 1/2 / — — —, y .' ' — -rexp(ixu)dx. A.14) 
V7 x i-oo x (cosh(7nc/4) - 1/2) PV 7 v ' 

Using the Euler-Maclaurin formula we can expand ( |A.10| ) and ( |A.11[ ), obtaining 
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and 

n / n / f f°° p(u — v)„... 1 , . . 1 p'(ii - A + ) 
<*(«) - aoo(n) = U -V^^* " 2iV^ " A+) + W^^T 



oo 



. , p(u + v) „ , . , 1 . . . 1 + A_) \ . . ^ „. 

+ -V^ ( ^ - ^ + A - } " w^cr j ' (A - 16) 

respectively. In the above equations A + and A_ are the largest and smallest root 
determined by the condition 

a n N {n)du = ^-(l+(3 ± (n)), (A.17) 
A± 2N 

where 

P±(n) = T^p- 

We should now consider separately the cases u > A + and u < — A_. In the case 
where u > A + (-u < — A_) the corrections of 0(1/ N 2 ) are calculated by neglecting 
the terms in the second (first), large bracket in (|A.16| ). Defining 

g(u)=p(u)/2iT, f ± (u) = a 00 (u + A±), x ± (u) = a n N {u + A±), (A.18) 

we can write ( |A,16| ) as 

This is precisely the standard form of the Wiener-Hopf equation (see, for exam- 
ple Morse and Feshbach ||16||). Its solution is obtained on defining the Fourier 
transforms 

roc I X ± (t) t><0 

xt(w) = / exp(iwt)x$(t)dt X ±(t) = { (A.20) 

J -°° t$0 

and the corresponding Fourier pairs g g, f f. Using the fact that 

(l-sMr^G+HG-H, (A.21) 

where 



r ( v ^r(l/2 - iw/A) expjiw ln(2)/4) 

G+H = r(5/e - »u,/8)r(i/6 - iw/s) = G - { - w) > (A - 22) 
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we can express x±( w )i after some algebraic manipulations as 

X+H = C±(w) + G + (w)(Q± + P ± (w)), (A.23) 

where 

C± (w) - 1 | ± 2 G + (i4/3)exp(-4/3A ± ) 

6 (W) " 2N + 12iV%- (A±) 4/3 -iw ' (A " 24) 

Equations (|A.19|) and the definitions (|A.18|) give us 
Q±(0) = —G+(iA/3) exp(-4/3A ± ) = — * gl - - - + (A.26) 

gf/2-4/3% % + 4/3 4/3/3±(n) 
fTjv(A±) " 12iVV^Q + 2N + 2iVG + (0)- (A - 27) 

Finally, using (|A.23|) , (|A.26| ) and ( |A.27| ) in ( |A.15| ) and approximating /(|u|) = 
v3exp(— 4/3|u|), u > 1, we obtain the first-order correction for the lowest energy 
= Ej^' in the sector with n dislocations and colour parity e, 

^ - e 3 ^ = ^5-(-l/6 + 2^) + o( l/iV 2 ), (A.28) 

where 

*< _ + <Lz£. ,A, 9 , 
In particular for the ground-state energy n = N, e = 1, we have 

§^^ = ^-^ + o(l/iV 2 ), (A.30) 



and comparing with ( |47| ) we obtain the value c = 1 for the conformal anomaly. If 
we now consider the gaps with respect to the ground-state with periodic boundary 
condition we obtain 

N N N 2 + l/ J ' 

Comparing this expression with (f48|) we obtain the conformal dimensions X^. 
These values are in perfect agreement with the operator content conjectured for 
the SCP and Baxter- Wu model, presented in section 4. 
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Figure and Table Captions 



Figure 1 - A triangular lattice with L = 3 rows and M=6 columns. The 
Baxter- Wu model is defined on the triangular lattice formed by the points 
o, □ and A, and the link variables {A} are defined on the hexagonal lattice 
formed by the points □ and o. The site-colouring problem (SCP) is defined 
on the hexagonal lattice formed by the points A and o. The open symbols o, 
□ and A denote the bordering sites related through (2) with the bulk ones 
(filled simbols). 

Figure 2 - Examples of configurations for the SCP with lattice size N = 4. 
The dislocations are represented by dotted lines. In (a) the configurations 
(m;X) = (1;2,3,4,4) and (m';X') = (3; 2, 2, 3, 4) are in the sector n = 
4 (re = 0). In (b) the configurations (m;X) = (1;2, 3,4) and (m';X') = 
(3; 2, 2, 2) belong to the sector n = 3 (re = 1). 

Table 1 - Ground-state energies per site for the SCP with lattice size N and 
boundary conditions re = 0, 1 and 2 given by (15). 

Table 2 - Conformal anomaly estimators c(N), as a function of the lattice 
size N, for the SCP and Baxter- Wu model. 

Table 3 - Scaling dimensions estimators Xj(N — n,l), as a function of the 
lattice size N, for some eigenenergies. These energies are the j th lowest 
energy obtained by solving fl50p-([52"D with values n, e and I. 

Table 4 - Energies for the Baxter- Wu model and SCP. The energies of the 
same line are identical. 

Table 5 - The mass-ratio estimators R,2{X,L) defined in (pop. 
Table 6 - The mass-ratio estimators R%{X,L) defined in fl60l) . 
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Table 1 



N 


K = 


K = 1 


K = 2 


6 


-1.3521881950 


-1.3396181681 


-1.3012024788 


10 


-1.3468393066 


-1.3423074850 


-1.3286264987 


50 


-1.3439405169 


-1.3437591407 


-1.3432148801 


100 


-1.3438498316 


-1.3438044868 


-1.3436684443 


150 


-1.3438330374 


-1.3438128841 


-1.3437524226 


200 


-1.3438271593 


-1.3438158231 


-1.3437818139 



Table 2 



N 


c(N) 


6 


0.996590995 


10 


0.998910268 


50 


0.999959561 


100 


0.999989915 


150 


0.999995519 


200 


0.999997480 
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Table 3 



N 


6 


10 


50 


100 


150 


200 


exact 


mo.o) 


0.12502803 


0.12501702 


0.12500083 


0.12500021 


0.12500009 


0.12500005 


0.125 


a* (1,0) 


0.24896741 


0.24959771 


0.24998323 


0.24999580 


0.24999813 


0.24999895 


0.25 


z+(0,l) 


0.50626226 


0.50215317 


0.50008406 


0.50002099 


0.50000933 


0.50000525 


0.5 


^(2,0) 


0.62613504 


0.62548322 


0.62502093 


0.62500524 


0.62500233 


0.62500131 


0.625 


xt(0,l) 


0.98648357 


0.99698967 


0.99992101 


0.99998057 


0.99999139 


0.99999516 


1 


z+(3,0) 


1.16335432 


1.13836199 


1.12552520 


1.12513123 


1.12505831 


1.12503280 


1.125 


zr(3,0) 


1.27059213 


1.25887795 


1.25038505 


1.25009649 


1.25004290 


1.25002413 


1.25 


4(0,2) 


1.53502476 


1.51231634 


1.50048724 


1.50012177 


1.50005411 


1.50003044 


1.5 


z+(0,l) 


1.82581639 


1.97119336 


1.99980095 


1.99995795 


1.99998195 


1.99998997 


2 



Table 4 



BW 


SCP 




E* UF {e = 1) 


E% w 


E^(e = -1) 




El CF {e = -1) 


E* w 


E^ F (e = 1) 




Ef F {e = 1) 


E BW 


E»<"(e = -1) 




E* UF {e = -1) 


E$ w 


E^ F (e = 1) 
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Table 5 



X 


7712 /mi 


2 


1.4049 


5 


1.7169 


6 


1.7273 


7 


1.7303 



Table 6 



X 


777.3/7711 


10 


2.02632 


20 


2.00538 


40 


2.00123 


40 


2.00053 
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Fig. 1 



{m\C'}=(3;2,2,2) 2 
{m,C}=(l;2,3,4) 



{m',C'}=(3;2,2,3,4) 
{m,C}=(l;2,3,4,4) 





